General Relativity and Gravitation manuscript No. 

(will be inserted by the editor) 



Phase transitions in geometrothermodynamics 

Hernando Quevedo • Alberto Sanchez ■ Sana 
Taj • Alejandro Vazquez 



Received: date / Accepted: date 

Abstract Using the formalism of geometrothermodynamics, we investigate the ge- 
ometric properties of the equilibrium manifold for diverse thermodynamic systems. 
Starting from Legendre invariant metrics of the phase manifold, we derive thermody- 
namic metrics for the equilibrium manifold whose curvature becomes singular at those 
points where phase transitions of first and second order occur. We conclude that the 
thermodynamic curvature of the equilibrium manifold, as defined in geometrothermo- 
dynamics, can be used as a measure of thermodynamic interaction in diverse systems 
with two and three thermodynamic degrees of freedom. 
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1 Introduction 

Although thermodynamics is based entirely upon empirical results which are satisfied 
under certain conditions in almost any macroscopic system, the geometric approach to 
thermodynamics has proved to be useful and illuminating. In very broad terms, one 
can say that the following three branches of geometry have found sound applications 
in equilibrium thermodynamics: analytic geometry, Riemannian geometry, and contact 
geometry. 

To illustrate the role of analytic geometry _1!;2_, let us consider a thermodynamic 
system specified by an equation of state of the form f(P, V, T) = 0. In the PVT— space, 
the equation of state defines a surface whose points represent the states in which the 
system can exist in thermodynamic equilibrium. Given a surface in a 3-dimensional 
space, its properties and structure are investigated in the framework of analytic ge- 
ometry. Probably, one of the most important contributions of analytic geometry to 
the understanding of thermodynamics is the identification of points of phase transi- 
tions with extremal points of the surface (or its projections on the coordinate planes) 
determined by the state equation f(P, V,T) = 0. More detailed descriptions of these 
contributions can be found in any good textbook on thermodynamics (see, for instance, 

EE])- 

Riemannian geometry was first introduced in statistical physics and thermodynam- 
ics by Rao [5], in 1945, by means of a metric whose components in local coordinates 
coincide with Fisher's information matrix. Rao's original work has been followed up 
and extended by a number of authors (see, e.g., [B] for a review). On the other hand, 
Riemannian geometry in the space of equilibrium states was introduced by Weinhold 
[7] and Ruppeiner 8,9, who defined metric structures as the Hessian of the internal 
energy and the entropy, respectively. 

Contact geometry was introduced by Hermann [TD] into the thermodynamic phase 
space in order to take into account the invariance of ordinary thermodynamics with 
respect to Legendre transformations, i.e., to the choice of thermodynamic potential. 
In order to incorporate Legendre invariance in Riemannian structures at the level of 
the phase space and the equilibrium space, the formalism of geometrothermodynamics 
(GTD) [TT] was recently proposed and applied to different systems. The main moti- 
vation for introducing the formalism of GTD was to formulate a geometric approach 
which takes into account the fact that in ordinary thermodynamics the description of 
a system does not depend on the choice of thermodynamic potential. One of the main 
goals of GTD has been to interpret in an invariant manner the curvature of the equi- 
librium space as a manifestation of the thermodynamic interaction. This would imply 
that an ideal gas and its generalizations with no mechanic interaction correspond to 
a Riemannian manifold with vanishing curvature. Moreover, in the case of interacting 
systems with non-trivial structure of phase transitions, the curvature should be non- 
vanishing and reproduce the behavior near the points where phase transitions occur. 
These intuitive statements represent concrete mathematical conditions for the metric 
structures of the spaces of phase and equilibrium. The main goal of the present work 
is to present metric structures which satisfy these conditions for systems with no ther- 
modynamic interaction as well as for systems characterized by interaction with first 
and second order phase transitions. 

This paper is organized as follows. In Section [2] we present the fundamentals of 
GTD. In Section[3]we analyze the case of systems with two and three degrees of freedom 
which are characterized by second order phase transitions. A metric is derived for 
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the equilibrium manifold which correctly reproduces the phase transitions in terms of 
curvature singularities. Section [4] is dedicated to the investigation of a thermodynamic 
metric which reproduces the first order phase transitions of realistic gases and fluids. 
Finally, Section [S] is devoted to discussions of our results and suggestions for further 
research. Throughout this paper we use units in which G — c = h — k B — 1. 



2 General framework 

The starting point of GTD is the (2n + l)-dimensional thermodynamic phase space 
T which can be coordinatized by the set of independent coordinates E a ,I a }, a — 
1, ...,n, where represents the thermodynamic potential, and E a and I a are the ex- 
tensive and intensive thermodynamic variables, respectively. The positive integer n 
indicates the number of thermodynamic degrees of freedom of the system. Let the fun- 
damental Gibbs 1-form be defined on T as = d<P—S a i > I a dE b , with 5 ab = diag(l, 1). 
The pair (T, 0) is called a contact manifold [101112] if T is differentiable and satisfies 
the condition A (d0) n ^ 0. Consider also on T a non-degenerate Riemannian metric 
G which is invariant with respect to Legendre transformations [13] 

(<l>,E a ,I a )^(3,E«,I"), (1) 

$ = i - 5 kl E k P , E i = —P, E j = E\ I 1 = E\ P = P , (2) 

where iUj is any disjoint decomposition of the set of indices {1, n}, and k, I = 1, i. 
In particular, for i = {1, n} and i = we obtain the total Legendre transformation 
and the identity, respectively. Legendre transformations are a special case of contact 
transformations which leave invariant the contact structure of T. Legendre invariance 
guarantees that the geometric properties of G do not depend on the thermodynamic 
potential used in its construction. The set (T, 0, G) defines a Riemannian contact 
manifold which is called the thermodynamic phase space (phase manifold), and repre- 
sents an auxiliary structure which is necessary to implement in a consistent manner 
the properties of ordinary thermodynamics. 

The space of thermodynamic equilibrium states (equilibrium manifold) is an n~ dim- 
ensional Riemannian submanifold £ C T induced by a smooth map tp : £ —¥ T , i.e. 
if : {E a ) i — y ($,E a ,I a ) with <P = $(E a ) such that (p*{0) = (p*{d<P - S ab I a dE b ) = 
holds, where tp* is the pullback of ip. This implies the relationships 

d$ = 5 ab I a dE b , ^ = S ab I b . (3) 

The first of these equations corresponds to the first law of thermodynamics, whereas 
the second one is usually known as the condition for thermodynamic equilibrium [3]. 
In the GTD formalism, the last equation also means that the intensive thermodynamic 
variables are dual to the extensive ones. Notice that the mapping <p as defined above 
implies that the equation <P = ${E a ) must be explicitly given. In standard thermody- 
namics this is known as the fundamental equation from which all the equations of state 
can be derived 3 , 4 . In this representation, the second law of thermodynamics is equiv- 
alent to the convexity condition on the thermodynamic potential d 2 &/dE a dE b > 0, if 
$ coincides with the internal energy of the system [TD] . The smooth map tp induces in a 
canonical manner a metric g (thermodynamic metric) in the equilibrium manifold £ by 
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means of g = <p*(G). Consequently, £ is a Riemannian manifold with a non-degenerate 
metric g. 

In general, the thermodynamic potential is a homogeneous function of its argu- 
ments, i.e., <&{\E a ) — \^$(\E a ) for constant parameters A and (3. Using the first law 
of thermodynamics, it can easily be shown that this homogeneity is equivalent to the 
relationships 

P$ = 5 ab I a E b , (1 - P)5 ab I a dE b + 5 ab E a dI b = 0, (4) 

which are known as Euler's identity and Gibbs-Duhem relation, respectively. 

From the above discussion we conclude that to describe a thermodynamic system in 
GTD it is necessary to specify a metric G for the phase manifold T and a fundamental 
equation $ — <P(E a ). These ingredients allows us to construct explicitly the pair (£,g) 
that constitutes the equilibrium manifold whose geometric properties should be related 
to the thermodynamic properties of the specific system. In the following sections we will 
present metrics G which generate equilibrium manifolds for different thermodynamic 
systems. 

3 Systems with second order phase transitions 

The non-degenerate metric 

G = (d® - S ab I a dE b f + (S ab E a I b )( Vcd dE c dI d ), Vab = diag(-l, 1, 1) (5) 

is invariant with respect to total Legendre transformations ([2| and, consequently, can 
be used to describe the geometric properties of the phase manifold T. The smooth 
map Lp : {E a } ^ {<P(E a ), E a , I a {E a )}, satisfying tp*{d<P - 5 ab I a dE b ) = 0, induces the 
thermodynamic metric 

'-•(°>- ('^al^^) < 6 > 

that can be explicitly calculated once the fundamental equation <P — ${E a ) is specified. 
Notice that by virtue of the equilibrium conditions Q and Euler's identity @, the 
conformal factor of the thermodynamic metric ([6| turns out to be proportional to 
the thermodynamic potential This means that in general we can assume that the 
conformal factor is different from zero. 

For concreteness we assume that E 1 — S is the entropy of a system with two 
thermodynamic degrees of freedom (n = 2). Then, the dual intensive variable I 1 = 
d<P/dS = T represents the temperature, and the thermodynamic metric (J6|) reduces to 

g = (ST + E 2 I 2 ) [-$ ss dS 2 + $22(dE 2 ) 2 ] , (7) 

where a subscript represents derivation with respect to the corresponding coordinate. 

To describe the phase transition structure of the system we use the heat capacity 
C = T(dT/dS) — $ s /$ ss , so that phase transitions of second order occur at those 
points where C diverges. For the calculation of concrete examples it is necessary to 
specify the fundamental equation $ = ^(S 1 , E 2 ). In particular, we are interested in ana- 
lyzing the thermodynamic properties of black hole configurations Q3] for which usually 
$ — M is the mass (energy) of the black hole and E 2 is an additional extensive variable 
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like electric charge, angular momentum, etc. In this case, the fundamental equation is 
essentially equivalent to the entropy-area relation S = S(M, E 2 ) = kA, where A is the 
area of the horizon and k is a constant that depends on the dimension of the space- 
time. However, it turns out that it is not always possible to express the fundamental 
equation in the form M = M(S,E 2 ) and the entropy-area relation S = S(M,E 2 ) 
must be used. In ordinary thermodynamics this corresponds to a change from the en- 
ergy representation to the entropy representation. GTD allows us to perform changes 
of representations in a simple manner. In fact, to obtain the entropy representation 
of the metric ([7]l we rewrite the first law of thermodynamics dM = TdS + I 2 dE 2 as 
dS = (l/T)dM - (I 2 /T)dE 2 and identify M and E 2 as the new extensive variables, 
whereas 1/T and —I 2 /T represent the corresponding new dual intensive variables. 
Replacing these new variables in the general expression (JS| , we obtain the metric 

g = (MS M + E 2 S 2 ) [S MM dM 2 + S 22 (dE 2 ) 2 ] , (8) 

which allows us to investigate the properties of the same thermodynamic system in the 
entropy representation. 

The above thermodynamic metrics (0 and (JS| have been used to investigate the 
black hole configurations listed in Table Q] In all the cases GTD is mathematically 
consistent and reproduces the thermodynamic behavior of the black holes. In fact, 
the scalar curvature of the equilibrium manifold (£,g) is in all cases different from 
zero, indicating the presence of non-trivial thermodynamic interaction. Furthermore, in 
Table[5]we present the heat capacity and the scalar curvature for each of the black holes 
contained in Table [T] It follows that at those points where the heat capacity diverges 
the scalar curvature of the equilibrium manifold becomes singular. Consequently, a 
second order phase transition is characterized by a curvature singularity. This shows 
that indeed the curvature of the equilibrium manifold can be considered as a measure 
of thermodynamic interaction. 

The fact that our results are invariant with respect to Legendre transformations ex- 
plains some contradictory results 21 , 22 23 that follow when the equilibrium manifold 
is equipped with metrics strongly associated to a particular thermodynamic potential. 

The 4-dimensional Kerr-Newman black hole represents a thermodynamic system 
with three degrees of freedom. Indeed, from the entropy-area relation one can derive 
the fundamental equation Q3] 

S = 7T [2M 2 -Q 2 + 2v/M 4 - M 2 Q 2 - J 2 ) , (9) 

which depends on the extensive variables M, J and Q. The second order phase tran- 
sitions are determined by the corresponding heat capacity: 

4TM 3 S 3 

r< — ^ al J (in) 

2M6-3M 4 Q 2 -6M 2 J 2 + Q 2 j2+2(Af 4 -i\/ 2 Q 2 - J 2 ) 3 / 2 ' K ' 

According to the general expression for the metric g, as given in Eq.©, the equi- 
librium manifold is 3-dimensional and the corresponding Legendre invariant metric 
reduces to 

MM 

9ab = ( MS m + QS Q + JSj) \ -s QQ -s QJ I I I I ! 

~Sq.j ~ "S 1 ././ 
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Black hole 



Fundamental equation 



RN 

RNAdS 
Kerr 

KAdS 
BTZ 
BTZCS 

BTZTF 
EGB 

EGBME 

EYMGB 
EMGB 



M . 



S = -n(M + M 2 - Q 2 ) 2 
£3-1 

167r y [i 2 + I 4S i + 



(D-2)(D-3)S 2 



5 = 2n{M 2 + VAf 4 - J 2 ) 



M - Ta^T? + "s^ - 



2 r2 \ D — 2 



M 



S 2 + 87T 2 fc J + f ^/(Z 2 - fc 2 )(S 2 + 167T 2 fcJ) 



S = 2V2ttI 



M - 



- § In 2^7^ 



Af - 



A/ 2 -£ 



M = 53 + 



M + 



V3 



M 



71 



M = S3 



as 3 



A/ + 2fiL 



s/3 



: lnS 



f53 -^53 



Table 1 Black holes with two degrees of freedom. Here we use the following notations: 
RN = Reissner-Nordstrom, RNAdS = Reissner-Nordstrom- Anti-de-Sitter, KAdS = Kerr-Anti- 
de-Sitter, BTZ = Bahados-Teitelboim-Zanelli, BTZCS = Bahados-Teitelboim-Zanelli-Chern- 
Simons, BTZTF = Banados-Teitelboim-Zanelli black hole with thermal fluctuations, EGB = 
Einstein-Gauss-Bonnet, EGBME = Einstein-Gauss-Bonnet black holes with modified entropy, 
EYMGB = Einstein- Yang-Mills-Gauss-Bonnet, EMGB = Einstein-Maxwell-Gauss-Bonnet, M 
= mass, 5=entropy, Q = charge, J = angular momentum, A = — l/l 2 = cosmological constant, 

J 1 ; 11 ' r l! ; ' ; - 2)-sphere, k = Chern-Simons 

(D - 2){D - 3)a. 



and oj£>^2 = 2-7T 2 j T[ — 2 1 ] is the volume of the unit (D 
coupling constant, a = Gauss-Bonnet coupling constant, a = 



Inserting here the expression for the entropy we obtain a rather cumbersome metric 
which cannot be written in a compact form. Nevertheless, the scalar curvature can be 



shown to have the form R 



KN 



KN it^KN 



where 



V KN = 2M 6 



3M 4 Q 2 



6M 2 J 2 +Q 2 J 2 +2(M 4 



M 2 Q 2 



r2x3/2 



(12) 



and M is a function which is always positive in the black- hole region M 4 > M 2 Q 2 + 
J 2 . Since the denominators of the heat capacity and the scalar curvature coincide, we 
conclude that there exist curvature singularities at those points where second order 
phase transitions occur. This result reinforces the interpretation of the curvature of 
the manifold £ as a measure of thermodynamic interaction. 

An additional test for the thermodynamic metrics and © consists in calculating 
the curvature of systems with no thermodynamic interaction or no phase transitions. 
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BH 



Heat Capacity 



Scalar curvature 



RN 



RNAdS 



Kerr 



KAdS 



BTZ 



BTZCS 



BTZTF 



EGB 



EGBME 



EYMGB 



EMGB 



ll^r 2 , (r + — r_) 
(r + -3r_) 



(D-2)S 



-^+(0-3) ^T3- J ~(-D-2) = 



~5~ 

"-1 - (0 - 3 )(^n °~ 2 ^ 2(2 °- 5) " 



27T r_[_ (r_j_ + r_ ) (r_i 



-) 



H — 6r_j_r_ — 3r_ 

(D-2)g[35 2 +2Q7r 2 J 2 -D(S 2 +47r 2 J 2 )](S 2 +4tt 2 J 2 ) 
35 4 + 24tt 2 J 2 S 2 + 2407r 4 J 4 -D(S 4 +487r 4 ,/ 4 ) 

47rr_j_ (t*5. — r ?_ ) 



47r(Z 2 -fc 2 )r 2 (r 2 -r 2 ) 
i[Zr + (r 2 + 3r 2 ) - fcr_ (r 2 +3r^ )] 



47rr_)_ (r^_ — r_ ) 

' r- 2 -L'Ar-' 2 ~ 



-35 | 3S 7~ q2 

3S3 -5Q 2 



2^1-4 / 9 2 -8i+(l+45)y / l-i| 2 - j 
S3 -Q 2 



S3 



35 



3S3 -AS^ 



5Q 2 -3S3 -AS 2 



(r + — 3r_ r++6r_ )(r_|_+3r_ )(r_^. — r_ ) 
,2^3 /„2 ,oJ -,2 f 



AT* 



-(£> 



-3)(^)' 



2(2D-5)ir 2 Q- ! 
(U-2)S2 



(3r^_+3r1_r_ +17r + r* +9ri ) (r + -r_ ) 3 
2ir 2 r 2 _(r + +r_) 4 (r 2 _-6r + r_-3r 2 ) 2 

.V 



' + 

rKAdS 



[3S 4 +24tt 2 J 2 S 2 +240ir 4 J 4 -D(S 4 +487r 4 ./ 4 )] 2 



3 ! 
3" 



2 (r 2 +3r 2 ) 
j^BTZCS 



- l|2 



! 2 [;r + (r 2 +3r 2 )-fcr_(r 2 +3r 2 )] 



47T 2 r 4 (r^_+3r 2 )3 



729 (3S 7 +10S V Q 2 -4s" Q 4 ) 
2 2S^(3S3-5Q 2 ) 2 (3S3+2Q 2 )3 



[l_4 / 32_ 8 i +( i +4 «5)] 2 (l- 4 |£) 

KfEYMGB 

2 \ 2 7 1 1 

Q2_ 1 S 3 j I Q 2 S3 +2Q 2 S3 In ; 

kj-EMGB 

4 W 4_ 

2Q 2 +3S3-AS 2 ) |5Q 2 -3S3-AS 2 



Table 2 GTD of black holes with two degrees of freedom. Here we use the following 
notations: r+ = radius of the exterior horizon, r_ = radius of the interior horizon, ft = Q/M, 
5 = a/M. The function AT represents in each case the numerator of the scalar curvature 
which is a well-behaved function at the points where the denominator vanishes. The results 
for charged rotating BTZ black holes are generalized in 1201 For more details see Refs. 15,16 



Table [3] contains the fundamental equations for some ordinary thermodynamic systems 
which we investigated from the point of view of GTD. 

Table [4] contains the results of our analysis. In the case of the ideal gas and its 
paramagnetic generalization, which have no thermodynamic interaction, the curvature 
vanishes and the equilibrium manifold becomes fiat. In fact, one can show the gen- 
eral result that any generalization of the ideal gas whose fundamental equation can 
be separate in its variables as S = Si(U) + S2CV) + Ss(E 3 ) + ... always generates a 
flat thermodynamic metric [M] . The 1-dimensional Ising model [55] generates a metric 
whose curvature is non-zero, indicating the presence of thermodynamic interaction, 



s 



System 


Fundamental equation 


IG 

PalG 

IIM 


S = JV(ln-£ + fln£) 
<J _ M []„ V , 3 lr , U 3 1rl 2M 2 \ 3 M 2 1 

ft - JV in w + ^ in ^ - 5 m [ i. - -j^r- 1-^ j 


F = — J — Tin ^coshf- + y/sinh 2 + e~ 4 T^j 



Table 3 Ordinary thermodynamic systems. We use the following notations: IG = ideal 
gas, PalG = paramagnetic ideal gas, IIM = 1— dimensional Ising model, U = energy, V = 
volume, N = total number of molecules, M = magnetization, H = magnetic field, T = 
Temperature, J = spin interaction parameter, F = Hclmholtz free energy. 



and regular everywhere, indicating the lack of second order phase transitions. Conse- 
quently, GTD reproduces at the level of the curvature the properties associated with 
the thermodynamic interaction between the particles of the mentioned thermodynamic 
systems. 



System 


Heat Capacity 


Curvature scalar 


IG 


3 TV 


R = 


PalG 




R = 


IIM 


1 f M(T,H) N(T,H) 2 \ 
\ j3[7+cosh f ] j2[/ +cosh Hf j 


R^O 



Table 4 GTD of some ordinary systems. Here / = I(T, H, J) is function with non-zero 
values. The Ricci scalar for the Ising model cannot be put in a compact form, but a numerical 
analysis shows that it is everywhere regular. 



4 Systems with first order phase transitions 

Consider the non-degenerate and Legendre invariant metric 

G = (d& - 8 ab I a dE b f + (5 ab E a I b )(S cd dE c dI d ), (13) 

for the phase manifold T. This metric is invariant with respect to total Legendre 
transformations @ and, consequently, can be used to describe the geometric properties 
of T. If we assume that the equilibrium manifold £ C T has as coordinates the extensive 
variables {E a }, the smooth embedding map £ T, satisfying tp*(d<P-S ab I a dE b ) = 0, 
can be used to generate in a canonical way the thermodynamic metric 

g = <p*(G) = ( h dE a dE b ) (14) 
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for the equilibrium manifold £ . In the case of a thermodynamic systems with only 
two degrees of freedom (n = 2), let us introduce the extensive variables <P = U = 
internal energy, E 1 = S = entropy, and E 2 = V = volume, together with the intensive 
variables T = temperature, and P = pressure. The first law of thermodynamics © 
reads dU = TdS — PdV and the explicit form of the thermodynamic metric is 

g = {ST- PV) (U ss dS 2 + 2U sv dSdV + U vv dV 2 ) (15) 

or, equivalently, in the entropy representation 

g = ^(U + VP) (S U[f dU 2 + 2S uv dUdV + S vv dV 2 ) . (16) 

Probably, the best-known thermodynamic system with a very rich structure of first 
order phase transitions is the van der Waals fluid [3J. For the sake of generality, we 
will use here the van der Waals fundamental equation together with the theorem of 
corresponding states in order to recast the fundamental equation in an invariant form 
applicable to all fluids 



S = ln(3V-l)+|ln(tf+|) , (17) 



where U = U/U c and V = V/V c , with U c = 4a/96, V c = 36, and a and b are the 
van der Waals constants. The heat capacity following from this fundamental equation 
turns out to be constant and, consequently, no second order phase transitions can 
occur. However, the critical points determined by the roots of the equation 

PV 3 - SV + 2 = 0, (18) 

correspond to first order phase transitions [3], where P = P/P c , with P c = a /27b 2 , is 
the reduced pressure. 

From the fundamental equation (|17l) it is then straightforward to compute the 
thermodynamic metric in the entropy representation: 

9 (wv 2 - UV -317 + 3) f_ 

g vdw = == — <WdU - 6V dUdV + 

4(3V - 1)(UV + S) 3 V I 



2(17 V + 3)(U V* - 6V 2 + 6V - 1) 
3(3F - l) 2 



dV 2 } . (19) 



The 2-dimensional equilibrium manifold turns out to be curved in general, indicating 
that the particles of the fluid interact thermodynamically. Furthermore, the scalar 
curvature of the above metric can be written in the form 

KfvdW 

R = ^ (20) 

t3 



( P V - 3V + 2 

where J\f vdw is a function of U, and V that is well-behaved at the points where the 
denominator vanishes. We see that the scalar curvature diverges at the critical points 
determined by Eq. ()18|) . Consequently, a first order phase transition can be interpreted 
geometrically as a curvature singularity. This is in accordance with our intuitive inter- 
pretation of thermodynamic curvature. 
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It is worth noticing that the metric (|14[) can also be used to describe the properties 
of ordinary systems. Indeed, a straightforward computation of the corresponding ther- 
modynamic metrics, using the fundamental equations of the ordinary systems listed in 
Table [31 leads to results equivalent to those reviewed in Table 2] 

5 Conclusions 

Geometrothermodynamics is a differential geometric formalism whose objective is to 
describe in an invariant manner the properties of thermodynamic systems in terms of 
geometric concepts. To a specific thermodynamic system, GTD assigns a particular 
equilibrium manifold £ whose geometric properties should be related to the thermo- 
dynamic properties of the system. In particular, we expect that the curvature of the 
equilibrium manifold is related to the thermodynamic interaction between the particles 
that constitute the system. 

In this work, we analyzed two different thermodynamic metrics which can be canon- 
ically derived from the non-degenerate Legendre invariant metric 

G = (d* - 5 ab I a dE b f + (5 ab E a I b )( Xcd dE c dI d ) (21) 

of the phase manifold T, where \ab ls an arbitrary constant diagonal tensor. If Xab — 
rjai), the resulting thermodynamic metric can be used to describe the properties of 
systems characterized by second order phase transitions. This was shown in particular 
for several black hole configurations in diverse theories and dimensions. In the Euclidean 
case, Xab — ^abi the corresponding thermodynamic metric was shown to correctly 
reproduce the structure of phase transitions of first order in the specific case of a 
fundamental equation which describes realistic gases and fluids. 

The formalism of GTD indicates that phase transitions occur at those points 
where the thermodynamic curvature is singular. The singularities represent critical 
points where the geometric description of GTD does not hold anymore and must give 
place to a more general approach. In ordinary thermodynamics the situation is simi- 
lar; near the points of phase transitions equilibrium thermodynamics is not valid and 
non-equilibrium models must be implemented. It would be interesting to study non- 
equilibrium thermodynamic models in the context of a generalized GTD. 

Our results show that the metric structure of the phase manifold T determines 
the type of systems that can be described by a specific thermodynamic metric: An 
Euclidean structure describes systems with first order phase transitions, whereas a 
pseudo-Euclidean structure describes systems with second order phase transitions. At 
the moment we do not have an explanation for this result. It is clear that the phase 
manifold contains information about thermodynamic systems; however, it will be nec- 
essary to further explore its geometric properties in order to understand where the 
thermodynamic information is encoded. In this context, the use of variational princi- 
ples could be of interest 
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